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We discuss the flavor conversion of neutrinos from core collapse supernovae that have oxygen- 
neon-magnesium (ONeMg) cores. Using the numerically calculated evolution of the star up to 650 
ms post bounce, we find that, for the normal mass hierarchy, the electron neutrino flux in a detector 
shows signatures of two typical features of an ONeMg-core supernova: a sharp step in the density 
profile at the base of the He shell and a faster shock wave propagation compared to iron core 
supernovae. Before the shock hits the density step {t ^ 150 ms), the survival probability of electron 
neutrinos above ~20 MeV of energy is about ~ 0.68, in contrast to values of ~ 0.32 or less for an 
iron core supernova. The passage of the shock through the step and its subsequent propagation 
cause a decrease of the survival probability and a decrease of the amplitude of oscillations in the 
Earth, reflecting the transition to a more adiabatic propagation inside the star. These changes affect 
the lower energy neutrinos first; they are faster and more sizable for larger ^13. They are unique 
of ONeMg-core supernovae, and give the possibility to test the speed of the shock wave. The time 
modulation of the Earth effect and its negative sign at the neutronization peak are the most robust 
signatures in a detector. 



PACS numbers: 97.60.Bw,14.60.Pq 

I. INTRODUCTION 

A core collapse supernova is mainly a neutrino phe- 
nomenon. The ~ 3 • 10^^ ergs of energy liberated in the 
collapse are radiated by neutrino emission from a thermal 
surface, the neutrinosphere. 

While the neutrino flux at the neutrinosphere is fairly 
independent of the star's properties, the neutrino flux we 
receive on Earth does depend on those. Indeed, neutrino 
oscillations are sensitive to the profile of the matter den- 
sity that the neutrinos encounter on their path. This 
makes them a very valuable tool to do a tomography of 
the star, in a way that depends on the neutrino mass 
spectrum and mixing matrix. An example is the possi- 
bility to use neutrinos to track the distortions in density 
caused by the propagation of the shock wave inside the 
star [T] . The sensitivity of the shock effects to the mixing 
9i3 would provide a powerful test of it |3 |3J HI |5J E] . 

An interesting implication of neutrino oscillations as 
tools for star tomography is the possibility to distinguish 
between different types of supernovae that differ in their 
density profile. In particular, neutrinos can distinguish 
[7] between a supernova with an iron core (Fe-core su- 
pernova) and one with an oxygen-neon-magnesium core 
(ONeMg-core supernova), for which the density distribu- 
tion and the shock propagation are completely different. 

Stars in the mass window between roughly eight and 
ten solar masses are expected to develop ONeMg cores, 
which may undergo gravitational collapse before Ne ig- 
nition due to rapid electron captures on ^^Mg and ^°Ne 
and may thus explode as so-called electron-capture super- 
novae (see for a recent study with a summary of previ- 
ous works). This distinguishes such stars from the more 
massive supernova progenitors, whose evolution through 
all stages of hydrostatic nuclear burning leads to the for- 



mation of an iron core at their center. ONeMg cores are 
bounded by an extremely steep density gradient (which 
we will call "density step" in the following), which dif- 
fers drastically from the much shallower density profiles 
around iron cores (see, e.g., figs. 1 of [21 110]) • 

The steep density decline allows the hydrodynamic 
bounce-shock to expand continuously, a fact that is fa- 
vorable for efficient neutrino heating in the post-shock 
layer and thus facilitates neutrino-driven explosions even 
in spherically symmetric simulations |1H I12j . Traveling 
down the density step, the outgoing shock also acceler- 
ates to much higher velocities than in the central regions 
of Fe-core supernovae. Despite the large difficulties in 
modeling their progenitor evolution, ONeMg core col- 
lapses have received quite some interest, in particular 
because of repeated speculations that they might be the 
long-sought site of the formation of r-process elements 
(e.g., [131 dU [121 Unj [T7]) and because of their possible 
link to the supernova of 1054 A.D., which gave birth to 
the Crab nebula (e.g., [HI HH |20l [211 122] ) • The number 
of similar stellar death events could account for a fair 
fraction of all supernovae. Poelarends et al. [H] estimate 
that about 4% of the stellar core collapses in the local 
universe might be of this kind, but the uncertainties in 
modeling the stellar evolution in the initial mass range 
between 6 and 12 Mq are significant and therefore the 
authors do not exclude a contribution of even 20% to all 
supernovae. 

To identify or exclude the oscillation signatures of an 
ONeMg-core supernova in an observed neutrino burst 
would be an important, unprecedented, test of stellar 
models and in particular of the existence of a density 
step in the interior of the star. It would be a necessary 
part of any data analysis aimed at reconstructing the 
originally produced neutrino fluxes (before oscillations) 
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and information on neutrino masses and mixings. 
But what are these signatures? 

Very recently Duan, Fuller, Carlson, and Qian have 
pointed out [7] that, due to the step in density, for 
ONeMg-core supernovae the neutronization peak (a peak 
in the Vf^ luminosity at about 10 ms post bounce) would 
not disappear regardless of the value of Q\-j,. In contrast, 
for just below the current bound disappearance due to 
v^-^v^^Vt conversion is expected for Fe-core supernovae 
[31[23J|21j. Duan et al. also pointed out how the conver- 
sion pattern is more complex for a ONeMg-core super- 
nova, because neutrino-neutrino scattering influences the 
resonant conversion driven by matter. The initial study 
in [7], as well as the analytical elaboration in [3S], were 
restricted to the neutrino flux at the neutronization peak 
and to relatively large 6'i3. Beyond these, the subject 
of neutrino conversion inside a ONeMg-core supernova 
remains to be explored. 

Here we give a in depth - even though far from compre- 
hensive - study of neutrino oscillations in ONeMg-core 
supernovae. We present the first discussion of shock wave 
effects for this type of supernova, using a numerical cal- 
culation of the density profile and of the neutrino fluxes, 
as they evolve over several hundreds of milliseconds. 

We describe a variety of conversion effects that could 
appear at different times depending on the value of ^la 
and on whether the detector is shielded by the Earth 
when the neutrino burst reaches it. We also address the 
question of what features could be observable and how. 
A first goal is to outline what combination of signatures 
one should look for to test models of ONeMg-core su- 
pernovae; another purpose is to describe the phenomena 
that should be taken into account when analyzing neu- 
trino data from a ONeMg-core supernova. 

The paper opens with generalities on ONeMg-core su- 
pernovae and details on our numerical model (secjiT]). In 
sec. |III| we give a discussion of conversion effects, followed 
by our results in sec. |IV[ A discussion on the implications 
of our findings closes the paper in sec. |V] 



II. OXYGEN-NEON-MAGNESIUM-CORE 
SUPERNOVAE: THE MODEL 

The interior structure of stars in the range ~8-10 Mq 
(Mq = 1.99 X lO'^^ kg is the mass of the Sun) is distinc- 
tively different from that of more massive supernova pro- 
genitors: their ONeMg core is surrounded by a thin car- 
bon and oxygen shell and an even thinner helium layer. 
The density gradient in this surface region of the core (at 
around 1000 Km) is extremely steep; the density can drop 
by nearly seven orders of magnitude within only 300 Km. 
This is visible from Fig. [T] where the line corresponding 
to time t = displays the electron number density pro- 
file of the innermost region of the 8.8 Mq progenitor star 
used as initial model for the core-collapse and supernova 
simulations this work is based on (baryon densities in the 
core and its surface layer can be obtained by multiplica- 



tion by a factor of 2) . The ONeMg core of this star was 
evolved to the onset of collapse by Nomoto [26l [27] and 
was more recently extended outside of the thin He shell 
by a hydrostatic hydrogen envelope pS] . 
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FIG. 1: Snapshots of the electron density profile at t = 
0, 50, 100,...., 700 ms (lines from top to bottom at their left 
end, except for the inverted curves for t = and t — 50 ms). 
The positions of the supernova shock for t > 200 ms coincide 
essentially with the lower right footpoints of the profiles. For 
t — 300 ms we also plot the eflective number density of neu- 
trinos (dashed curve, see Eq. Q), which is responsible for 
the effects of neutrino-neutrino forward scattering. The two 
horizontal lines represent the densities corresponding to the 
two MSW resonances for a neutrino of 20 MeV energy. 

The updated progenitor model was followed through 
collapse and explosion with the neutrino-hydrodynamics 
code Vertex [29l [30]. Details of the employed input 
physics and information about the supernova dynamics 
can be found in the works by Kitaura et al. [12, and 
Janka et al. [ni[TD]. Figure [l] provides an overview of the 
evolution of the electron density (rie) profile in steps of 
50 ms until 700 ms after the start of a spherically sym- 
metric (ID) simulation. The core bounce happens at 
t = 53.6 ms and the explosion sets in at about 130 ms (to 
be recognized from positive velocities developing in the 
post-shock layer). 

For t = 100 ms the step in density induced by the shock 
is visible at r ~ 130 Km. The shock front advances 
outward rapidly, reaching r ~ 500 Km at i < 150 ms 
and meeting the base of the He shell (at r ~ 1100 Km 
and ^ 4: X 10^^ cm^^) less than 10ms later. At this 
point, the shock accelerates even more and travels essen- 
tially with the speed of light until it is gradually decel- 
erated again by running into the flatter density profile of 
the hydrogen envelope and sweeping up matter from the 
progenitor star. At the end of the computed evolution, 
the shock has reached a radius of more than 2x10'"' Km 
(Fig.0. 

As a consequence of the differential collapse of the core, 
starting with the innermost regions and only gradually 
encompassing layers farther away from the center, the 
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density step is significantly flatter at the time when the 
shock hits the ONeMg-core surface and passes through 
the base of the He shell. For the first roughly 100 ms af- 
ter the shock passage through the step, the density pro- 
file therefore appears essentially smooth. The ejecta in 
the region of the previous density step then expand with 
nearly uniform and constant velocity (~2-3x 10"* cms~^) 
and the width of this layer (1000-2000 km) is preserved. 
At later times, however, the profile appears to steepen 
again because of the dramatically increasing radial scale 
and the resulting relative decrease of the ratio of width to 
radius of the shell. The relics of the original density step 
of the pre-coUapse core can therefore still be recognized 
as an outward moving step-like structure on the density 
profile of the shock-accelerated and expanding ejecta of 
the supernova explosion. 
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FIG. 2: Neutrino luminosities (top) and average energies (bot- 
tom) before oscillations. The Vx (vx =^'^1 ^t, J^/j, I'r) luminos- 
ity is per species, so the total luminosity in the non-electron 
flavors is four times the one plotted here. The average energy 
is defined as the ratio of energy flux to number flux. 

The time evolution of the neutrino luminosities and 
of the mean spectral energies (measured at a radius of 
r = 400 Km for an observer in the laboratory frame), 
without including flavor conversion, is shown in Fig. [2] 
The displayed data are from a two-dimensional (2D) sim- 
ulation (assuming axisymmetry around the polar grid 
axis) of an ONeMg-core supernova, which was followed 



until roughly 200 ms after core bounce. The explosion dy- 
namics of the ID and 2D models is very similar and the 
density profiles as shown in Fig. [T] basically also apply for 
the 2D model after averaging over latitudes (thus wiping 
out inhomogeneities associated with convcctive overturn 
in the neutrino- heating layer around the neutron star). 
In contrast to ID results, the neutrino luminosities and 
mean energies from the 2D model include the modifica- 
tions caused to the neutrino emission by convective en- 
ergy transport inside the nascent neutron star. This plays 
a role later than a few ten milliseconds after core bounce 
and leads to slightly (about 10%) enhanced electron neu- 
trino and in particular muon and tau neutrino luminosi- 
ties, a slightly reduced electron antineutrino luminosity, 
and somewhat lowered ('^O.S-lMeV) mean energies of 
all radiated neutrinos ^31] . 

Figure |2] shows that the luminosities, flavor composi- 
tion, and mean spectral energies of the neutrinos reveal 
the well known behavior: the peak in associated with 
the neutronization burst, the slower rise of the Pg and 
iyx =i^/ij i^Ti ^•^) luminosities, and the hierarchy 
of the !/e, Ve and average energies. The main differ- 
ence compared to Fe-core supernovae is the absence of an 
extended phase of post-bounce accretion of the forming 
neutron star and therefore the lack of the corresponding 
accretion plateau in the neutrino luminosities after the 
Vf. shock-breakout burst. 

The 2D simulation was only carried to about 200 ms 
after bounce. For our analysis of neutrino oscillation ef- 
fects at later times until the end of the ID simulation at 
~ 650 ms post bounce, we will take the neutrino luminosi- 
ties and mean energies to be roughly constant. This is 
a reasonably good assumption because of the mentioned 
absence of accretion and the only slow, quasi-steady evo- 
lution of the settling neutron star. Moreover, our discus- 
sion will mostly concentrate on matter-of-principle effects 
and the details of the neutrino emission properties will 
not be relevant. 

The fluxes in Fig. [2] will be used as initial conditions 
for the calculation of neutrino oscillations, making the 
assumption that the neutrino luminosities and mean en- 
ergies are constant outside of the neutrinosphere. 

In our calculations, we model the neutrino energy spec- 
tra following Keil, Raffclt and Janka [32] : 

dN ^ {l + aY+-L (E\ " 

dE V{\ + a)Eo'\Ej ' 

where Eq is the average energy and L the luminosity. We 
take a = 2.3, for which the spectrum is close to a Fermi- 
Dirac distribution. The power law in Eq. ([ij is adequate 
in the fact that it reproduces well the numerical results 
of sophisticated neutrino transport (though by choosing 
a constant value for a we neglect the detailed evolution 
of the spectral shape provided by the simulation results) 
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III. FLAVOR CONVERSION IN THE STAR 
AND IN THE EARTH: GENERALITIES 

Neutrino conversion in supernovae is particularly com- 
plex due to the interplay of several effects. Here we re- 
view them briefly, in the measure needed to highlight 
the phenomena that are distinctive of a ONeMg-core su- 
pernova. We refer to the literature for more complete 
reviews (e.g., [H IMJ IM])- 

The flavor evolution of neutrinos inside the star and 
in the Earth is described by a Hamiltonian that contains 
three terms. The first is the kinetic term, which depends 
on three mixing angles 9ij, a CP- violating phase Sep, and 
two independent mass squared splittings, Am|j, where 
the indices i,j = 1, 2, 3 refer to the three mass eigenstates 
J^i, i^2) ^'s- The three characteristic oscillation frequencies 
associated with the kinetic term are: 



Ar 



2E 



(2) 



with E the neutrino energy. We use the standard 
parametrization of the mixing matrix, with the values 



sm 
10- 



Cl2 — 

3 eV2, 



0.32, 



and sin^ 9i3< 2 



8 • 10-'^ eV2, 

10~^ (see e.g., the summary 



lAm^il 



in [35] and references therein). Notice that 6*13 has not 
been measured yet [351 ISZ] , and that Ato§^ is known only 
in absolute value. The case of positive (negative) Am'^i is 
called normal (inverted) mass hierarchy. The phase 6c p 
can be ignored here, as it does not influence the relevant 
probabilities for equal i/^ and fluxes [3S1 ISH] ■ 

The second and third terms relevant for neutrino con- 
version are due to neutrino forward scattering (refrac- 
tion) on electrons and on neutrinos. Their characteristic 
frequencies are 



= V2G 



pne 



= V2GFnl- 



eff 



(3) 



where Ue is the number density of electrons along the tra- 
jectory of the neutrinos, and nl/-^ is the effective number 
density of neutrinos. It is deflned as (see e.g., [M]): 




(4) 



where r is the radial coordinate, R is the radius of the 
neutrinosphere, and N^, is the total number density of 
neutrinos and antineutrinos of all flavors at the neutri- 
nosphere: = J2a=e fi ri^a^^a)- The self-interactiou 
term from neutrino-neutrino scattering generates non- 
linear effects that only recently have been studied in de- 
tail in the context of supernovae [71113 [Ml HOI HD 113 
SHI mi Hg [4611471 SHI [49]. 

The character of the neutrino conversion is roughly 
determined by the relative size of the different frequen- 
cies ujji, uje, and lo^. Each of the two refraction terms is 
relevant if its frequency is comparable or larger than the 
kinetic (vacuum) ones. The more terms contribute simul- 
taneously, the higher is the complexity of the conversion 



pattern. Fig. [T| allows to compare the strengths of the 
vacuum and refraction terms. Together with the elec- 
tron density proflle at different times, the figure shows 
the values of n,, for which the conditions LO21 — uJe and 
Lo^i — LOf. are realized for E — 2^ MeV. For small mixing 
angles (as it is the case for 6*13), these well approximate 
the condition for the MSW resonance [H [SS [SI] [5^ . We 
also plot the effective neutrino number density, nlJ^ (r) 
for t = 300 ms and i? = 60 Km. This quantity is smaller 
(larger) at later (earlier) times, as follows from fig. |2] 

From fig. [1] we identify two different scenarios: the 
stage before the shock reaches the He shell, t ^ 150 ms, 
and the later times, t 150 — 200 ms. These phases 
will be called pre-shock and post-shock respectively 
from here on. In what follows we discuss them for the 
propagation of neutrinos with the normal mass hierarchy, 
which is the most transparent and sufficient to illustrate 
the main effects. 



A. The late stage: post-shock 

Let us first consider the late, post-shock phase, for 
t ^ 300 ms. For this, in the most internal part of 
the star (r ~ 100 Km) both the neutrino-electron and 
the neutrino-neutrino scattering terms exceed the vac- 
uum one and are comparable to each other. In general, 
they both affect oscillations substantially. As the neu- 
trinos propagate out to larger radii, the influence self- 
interaction is strong over a typical distance of 10^ — 10^ 
Km. At r ~ 3000 Km, the neutrino-neutrino term, 
ujiy, falls below the vacuum frequency LO21, and thus be- 
comes negligible, marking the end of self-interaction ef- 
fects. Since the electron number density decreases more 
slowly with the radius, at this distance the neutrino- 
electron term is still large and dominates over the vac- 
uum one by at least one order of magnitude. Only at 
radius r ^ 4 — 5 • 10"^ Km we have 0^31 ~ Wg, which drives 
the inner MSW resonance. Another resonance, driven 
by UJ21 ^ ^e, follows at even larger radius, r ^ 8 • 10^ 
Km, beyond which all matter effects are negligible. The 
separation between the distance of decoupling of the self- 
interaction effects and the position of the MSW reso- 
nances, as well as the separation between the two res- 
onances themselves, increases with time. Thanks to 
this separation, the phase of collective effects, driven by 
neutrino-neutrino scattering, and the MSW resonances 
are decoupled. 

A scenario with these characteristics has been studied, 
for a Fe-core supernova, by Fogli et al. ^34]- They found 
that, beyond the region where self-interaction effects end, 
the neutrino spectra are unchanged for the normal mass 
hierarchy. Such conclusion has little dependence on ^13 
and on the speciflc electron density profile, as long as 
the condition We ^ W31 is satisfied over the distance of 
effectiveness neutrino-neutrino scattering ^4] (see also 
[121 H5] for further discussions) . It is supported by ana- 
lytical arguments based on the minimization of potential 
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energy of a spin in a magnetic field (e.g., [Ml Hll ITT] ) . 
Only with the violation of the matter-dominance condi- 
tion (which is not our case) a swap of the v^, and / Vt 
spectra might occur below a critical energy Ec ^ 10 
MeV [3B] . Even in the implausible case that a swap with 
Ec ^ 10 MeV occurs in our scenario, such swap would 
depend only indirectly on the electron density, with no 
strong sensitivity to the density structure of a ONeMg- 
core supernova. 

For all these reasons, we consider it adequate here to 
consider unchanged neutrino spectra immediately after 
(in radial distance) the end of self interaction effects, and 
focus of the physics of the MSW-driven transformation, 
which, thanks to its direct and resonant dependence on 
Ue, is a much more sensitive probe of the matter distribu- 
tion inside the star. As will be seen, the most interesting 
effects of the MSW resonances are at high energy {E > 30 
MeV) , and thus are not influenced by a possible swap at 
low energy. 

After the MSW-driven conversion, the flux of electron 
neutrinos at Earth, F^,, can be written in terms of the 
fluxes without oscillations, F^ and F° {x = r), as: 



(5) 



where p denotes the Vg, survival probability. This proba- 
bility is given by 



P~Ph [cos^ 0i2Pl + sin^ ^12(1 - Pl)] 



(6) 



Here Ph {Pl) is the probability of transition between 
the eigenstates 1^3 and 1^2 {1^2 and vi) of the Hamiltonian 
in the higher density (lower density) MSW resonance. 
Ph depends on ^13 and on the derivative of the electron 
number density, calculated at the point Vres where W31 = 
uje (see e.g., [53]): 
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(7) 



Pl obeys the same expression with the replacements 
wai 0-121, 6*13 012 and Arri^i Am|]^. We can see 
that the survival probability varies between p = and 
p = cos^ 012 (numerically, between p = and p = 0.68), 
for Pff and Pl varying in the allowed range from (adi- 
abatic propagation) to 1 (completely non-adiabatic prop- 
agation) . The derivative of the density profile in Eq. Q 
encodes the sensitivity of the MSW effect to the matter 
distribution: the amount of flavor conversion depends 
on how rapidly the electron density changes along the 
neutrino trajectory and on the size of the mixing an- 
gle. Steeper profiles and smaller mixings correspond to 
stronger violation of adiabaticity. 

The post-shock phase can be divided in an early part, 
t <, 650 ms, and a late one, t ^ 650 ms. The late part 
begins when, accelerated by the shock, the He shell has 



traveled outwards and expanded to lower density, to a de- 
gree that the density step associated with it does not ex- 
tend to the resonance density anymore, and thus becomes 
irrelevant for neutrino conversion. In coincidence with 
this transition to a shallower density profile, we expect 
the MSW conversion to become more adiabatic. Such 
early transition in the direction from less to more adia- 
batic conversion is unique of a ONeMg-core supernova. 



B. The early stage: pre- shock 

Let us now discuss the more complicated pre-shock 
phase. Here, all the three frequencies are comparable 
at r ~ 10^ Km. Therefore the physics of the MSW 
resonances can not be decoupled from that of neutrino- 
neutrino scattering as in the post-shock phase. Moreover, 
the two resonances are spatially closer, to the point that 
the factorization of their effects may not be valid. 

In ref. [7] the Ve survival probability for this case was 
calculated numerically, for ^13= 0.1. It appears that, for 
the normal hierarchy, the probability can be described by 
transitions between the eigenstates of the Hamiltonian, 
in a way that resembles the simpler case of two spatially 
separated resonances without effects of neutrino- neutrino 
refraction. An important feature is the presence of sharp 
swaps in the energy spectrum, instead of the smoother 
transitions between adiabatic and non-adiabatic regimes 
that are expected in the pure MSW case as the neutrino 
energy increases. 

In the limit of slowly varying neutrino density, the 
probability found in 7 has an analytical interpretation 
[5S] : it is the result of collective MSW transitions - where 
the collective behavior is due to neutrino-neutrino scat- 
tering - followed, at lower density, by spectral swaps anal- 
ogous to those already observed for other density profiles 
[40l |46] . The survival probability can be described 
with the same expression, Eq. Q, using effective values 
for Ph and Pl [6S]. In first approximation, these behave 
as step functions, whose critical energies are determined 
by conservation laws ^25^ i49j . Specifically, the numerical 
result in [7] is reproduced by: 



Ph = 



Pl 



) 1 for £: > 12 MeV 

1 otherwise 

'1 for£'>15MeV 

otherwise 



(8) 



so that we have p ~ below 12 MeV and p ~ 0.68 
at i? ^ 15 MeV, with p ~ 0.32 as intermediate value 
between the two. 

While the dependence of the results on 6*13 was not 
shown in [7], it was checked by the same authors [SI] that 
the behavior p = cos^ 812 — 0.68 at high energy (above a 
critical energy Ec ^ 15 MeV) is a generic feature even 
for smaller 6*13, reflecting the fact that the produced z^e's 
always emerge the star almost completely in the lightest 
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mass state, vi. This can also be understood considering 
that in the hmit ^13^ 0, the third mass state is a pure 
mixture of i^^ and Vt, and therefore it decouples from 
the evolution of z/g. This implies that the I'e survival 
probability should depend only on the pair of states 1^2 
and z^i . These have a content ranging from sin'^ 612 = 
0.32 to cos^ 612 = 0.68. A 1^2 ^ i^i swap would produce a 
change in the survival probability from the lower to the 
higher of these two values. 

Here we use the numerical result from ref. ^ for the 
survival probability. For ^13= 0.1, it applies well to our 
case because all the parameters match: our pre-shock 
density profile matches the one in [7] and the neutrino 
spectrum at the neutronization peak (t ~ 60 ms) has a 
very similar average energy (12-13 MeV against the 11 
MeV in [7]) [SS]. For smaller ^13, we consider the same 
result to be still applicable to the high energy range, in 
the light of the reasoning above. As it will appear shortly, 
the high energy range is where we expect the most robust 
signatures of conversion, namely the Earth regeneration. 

Notice that we do not address the neutrino conver- 
sion at the instants of time that immediately precede 
and follow the shock passage through the density step, 
t ~ 150 — 300 ms. This is because in this interval the con- 
version pattern is highly complicated: the matter near 
the resonance density undergoes a quick acceleration to 
relativistic velocity with the passage of the shock, thus 
making relativistic corrections necessary [SS]. Instead, 
for the same density, velocities are lower than 10% of the 
speed of light at later times. In addition to relativistic 
corrections, one would need to include all the oscillation 
terms (kinetic, electron scattering and self-interaction) 
simultaneously, as in the pre-shock phase. The combina- 
tion of these effects would require a dedicated study that 
is postponed for the time being. 

C. Oscillations in the Earth 

In addition to conversion in the star, neutrinos from a 
supernova undergo oscillations inside the Earth. These 
oscillations are observable by a detector that is shielded 
by the Earth when the neutrino burst reaches it. They 
are driven by neutrino-electron scattering. The phase 
of oscillation depends on the solar parameters, 612 and 
Am2i, while the amplitude reflects the pattern of con- 
version in the star. The difference between the Ve flux 
in a detector with and without Earth shielding has the 
expression (valid for normal hierarchy) I56j: 

- Fe = (F," - F,")P^(1 - 2PL)freg ■ (9) 

Here freg is the regeneration factor. Up to terms propor- 
tional to sin^ ^13, freg is given by: 

/,eg=P®(i^2->i^e)-sin2 012 , (10) 

where P®{y2 -^Vf.) is the probability that a state enter- 
ing the Earth as V2 is detected as in the detector. In 



the absence of Earth shielding, P®{u2 -^^e) = srn^ 9x2 
and freg = 0. The regeneration factor is an oscillatory 
function, and freg ~ with good approximation at all en- 
ergies and zenith angles. Notice that Eq. (|9| exhibits an 
elegant factorization of the several steps of the neutrino 
propagation from production to detection: the first fac- 
tor describes the originally produced fluxes, the second 
the higher density resonance, the third the low density 
resonance and the fourth the effect of the Earth. This 
factorization makes the Earth effect particularly trans- 
parent, and the oscillatory character makes it an unam- 
biguous signature of neutrino flavor conversion, that can 
not be confused with astrophysical effects. 

Here we calculate the oscillations in the Earth using 
Eq. ([9| with and as described above for the two 
regimes (pre and post shock), and with freg from the 
accurate numerical calculation in [56] . 

IV. CONVERSION OF NEUTRINOS FROM AN 
ONEMG-CORE SUPERNOVA: RESULTS 

A. Probabilities 

From the previous sections, one can see two oscillation 
effects that are distinctive of a ONeMg-core supernova. 
One of them is the influence on the adiabaticity of con- 
version of the sharp gradient in density that marks the 
oxygen-helium transition in the star. One expects strong 
adabaticity breaking even for relatively large Oi^, in con- 
trast to Fe-core supernovae. The second effect is the time 
variation of the oscillation probabilities within the flrst 
hundreds of milliseconds of the burst, associated with 
the variation of the density proflle with the shock pas- 
sage. In Fe-core supernovae shock effects appear at much 
later times. Below we illustrate these effects in detail and 
compare them to the results for Fe-core supernovae from 
the hterature (e.g., [531 [321 133) ■ In our calculations we 
have used the results of Duan et al. for the pre-shock 
phase, and the density profiles in fig. [l] together with 
the analytics in sec. |III| for the post-shock times. 

Let us begin with the jumping probability Pp. In fig. 
|3]Pff is plotted as a function of sin^ 6*13 for F = 20 MeV 
and for three representative times, one in the pre-shock 
phase (t = 60 ms [70]), one in the early post-shock phase 
(t = 450 ms) and one in the late post-shock stage {t = 700 
ms) . The same function for a Fe-core supernova (at i ^ 3 
s, from ref. [5]) is also shown for comparison. 

We see essentially three regimes in ^13: 

1. large: sin^ ^13 ^ 3 • 10"'^. The character of the 
conversion changes fast (within ~ 200 — 300 ms or 
so), turning from non-adiabatic in the pre-shock 
phase, to completely adiabatic afterward. 

2. intermediate: sin^ 6*13 ~ 3 • 10"^ - 3 • 10"^ The 
character of conversion changes more slowly: it is 
completely non-adiabatic in the pre-shock phase, 
and after ^ 300 ms it is still at least partially 
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FIG. 3: The jumping probability Ph for t = 60, 450, 700 ms 
(solid curves, from upper to lower) as a function of sin^ S13 
for energy E = 2Q MeV. The dashed line shows the same 
probability for a Fe supernova with the parameters in ref. [3]. 



non-adiabatic. Only in the late post-shock phase 
(t ^ 700 ms) it may become completely adiabatic 
if sin^ ~ 3 • 10"'' — 3 • 10"'^. In this range of mix- 
ing the three time intervals (pre-shock, post-shock 
early and post-shock late) can be distinguished. 



3. small: sin 6*13 ^3-10 



The conversion remains 
completely non-adiabatic at all times, with only a 
minor change in the transition from the early to the 
late post-shock phases. 

Notice that only in the late post-shock stage the de- 
pendence of Ph on ^13 resembles the one of a Fe-core su- 
pernova, for which there is no time variation until much 
later times. 

The jumping probability in the low density resonance, 
Pl, is found to be zero at all post-shock times, and is set 
to 1 before the shock for E ^ 12 MeV, as discussed in 
III For a Fe-core supernova Pl = at all times. 
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FIG. 4: The electron neutrino survival probability for t = 
60, 450, 700 ms (solid, short dashed, and long dashed, respec- 
tively) and sin'^ 613 = 6 ■ 10"^. p = for a Fe-core supernova 
over the same time interval. 



The behavior of P/,, combined with the results in fig. [3] 
for Ph, explains what we see in fig. [4j showing the energy 



dependence of the survival probability p for sin ^13 = 
6- lO"** and no Earth shielding. In the transition from the 
pre- to the post-shock phases the probability decreases 
from ^ 0.68 to ~ 0.32 (asymptotic values at high energy) 
mostly as a consequence of the change in P^. At later 
times, p goes to zero, beginning with the lower energy 
part of the neutrino spectrum. This appears in fig. [4] for 
t — 700 ms below E ^ 22 MeV we have p = 0, because for 
these energies the high density resonance is realized in the 
shallower part of the density profile, more internal with 
respect to the base of the He shell (regrouped after the 
shock passage). Instead, at higher energy the resonance 
is realized at lower density, on the density step. The 
step is moving outwards and becoming less dense, so the 
transition to p ~ moves to larger energy. 

For small 613, Ph — 1 at all times, and so the only 
change in p is due to the transition of P^ : the asymptotic 
value of p goes from ^ 0.68 to ^ 0.32 and retains this 
value at all later times. 

For large 6'i3, the probability p changes from p ~ 68 to 
zero immediately after the shock passage, when Ph goes 
to zero. 
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FIG. 5: The Ve survival probability as a function of time for 
E = 20 MeV and sin^ 6I13 = 0.01, 6 • 10"", 10"^ (solid, short 
dashed, and long dashed, respectively). The step-like struc- 
ture reflects the fact that the oscillation effects were evaluated 
for density profiles in steps of 50 ms (t = 0, 50, 100, .. ms, mid- 
dle points of the steps, see profiles of fig. [T]). For t ~ 150 — 300 
ms the figure is only schematic: it does not capture the full 
complexity of the conversion in accelerating matter when the 
shock passes through the density step (see sec. [Ill B[ ). 

The features seen in fig. [4] appear also in fig. [5] where 
the time dependence of p is shown for fixed energy and 
different values of ^13. The figure is a good quantitative 
description for the pre-shock and late post-shock phases 
(t 300 ms), while for the intermediate times it has 
only a schematic character. This is because it does not 
include a number of effects that are relevant when the 
matter at the base of the He shell is strongly accelerated 
by the shock (see sec. IIIB I. The figure evidences clearly 
the general decline of the survival probability with time, 
up to minor deviations from the monotonic trend due 
to the minute details of the density profile at the MSW 
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resonance. It is also apparent that this decrease is faster 
for larger ^13, as already noticed above. 

The conversion pattern found for a ONeMg-core su- 
pernova is in contrast with that for a Fe-core supernova 
over the same time interval: in that case, the survival 
probability has no time dependence and ranges from 
p — svo? O12 — 0.32 to p — Q depending on 6*13 (fig. [s]). 

In the absence of Earth shielding, the observation of 
the oscillation effects could be challenging. Indeed, the 
measurement of a probability would be complicated by 
uncertainties in the original neutrino fluxes, and the 
shock-induced modulations of the probabilities could be 
masked by the natural time evolution of the neutrino 
spectra and luminosities. The best signature to look at 
would probably be the fate of the peak in the Ve luminos- 
ity at i ~ 60 ms: for large ^13 it disappears for a Fe-core 
supernova, but it survives for an ONeMg-core one [7]. A 
conclusion on this would require knowing ^13, however. 
Looking for the progressive vanishing of p from lower to 
higher energy (fig. |4| might also be promising. 

As it has already been pointed out for a Fe-core su- 
pernova [571 [5^ , the Earth shielding can be a great 
advantage, because oscillations in the Earth have an un- 
ambiguous signature: oscillatory distortions in the ob- 
served neutrino energy spectrum. The very presence of 
the effect allows one to conclude on neutrino masses and 
mixing, even in the presence of uncertainties on the orig- 
inal neutrino fluxes. 
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FIG. 6: The relative Earth matter effect in the neutrino chan- 
nel, at 60 degrees nadir angle, for t = 60,450, 700 ms (solid-, 
short-dashed and long-dashed, respectively) and sin^ On = 
6 • 10~*. For a Fe supernova and the same value of ^13, the 
Earth effect is zero. 

Fig. ^shows the relative Earth effect, defined as {F^ — 
Fe)/i^e (see Eq. ([9])), for nadir angle of 60 degrees, an 
intermediate value oiOi^ (sin^ 6*13= 6-10~^) and the same 
instants of time as in the previous plots. The amplitude 
of the effect reflects the time evolution of Ph and Pl ■ In 
the pre-shock regime the amplitude is maximal, reaching 
- 25% effect at £: ~ 56 MeV. The effect is negative 
because, with Pl ^ 1, the term (1 — 2Pi) is negative and 



the flux contribution is positive, — > 0, reflecting 
the fact that Ve dominates the flux at this stage. 

In the early post-shock phase the effect is reduced in 
amplitude at low energy, where the conversion is partially 
adiabatic in the high density resonance. Notice that, with 
respect to the pre-shock phase, above £^ ~ 15 MeV there 
is a double change of sign in Eq. ([o]): one is caused 
by Pl changing from 1 to 0, and the other due to the 
transition from dominated to i^aj-dominated (at high 
energy) flux, meaning a change of sign in the term F^ — 
F^. The net Earth effect is again negative. Notice that 
the double change of sign is a unique feature of a ONeMg- 
core supernova: for any value of ^13 at t = 60 ms the 
relative Earth effect would be positive in a Fe supernova. 

We stress that the two sign flips may occur at slightly 
different times, thus giving rise to an interesting sequence 
of sign changes in the Earth effect. It might also be possi- 
ble to see the amplitude of the oscillations become smaller 
and then increase again if the transition of Pl from 1 to 
is gradual and Pl ~ 1/2 for a significant period of time 
and/or if the quantity F^ — F^ reaches a minimum in 
absolute value. To study if and how these time modula- 
tions happen is beyond the scope of the present paper: 
it would require information on the time evolution of the 
density profile with smaller time steps, and to use the 
time evolved profile for a detailed numerical study of the 
neutrino conversion in a regime where all the three oscil- 
lation terms are relevant. 

In the late post-shock phase the Earth effect is fur- 
ther reduced in amplitude and ultimately goes to zero 
together with Ph, starting with the lower energy part of 
the neutrino spectrum and including higher energies as 
time passes. In Fig. [6]we see the same feature at = 22 
MeV as in fig. |4] marking the transition to completely 
adiabatic conversion in the high-density resonance at late 
post-shock times. 
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FIG. 7: Same as fig. |6]for sin^ S13 = 10~^. The conversion 
is non-adiabatic at all times, and so the oscillations in the 
Earth never vanish. They differ from the case of a Fe-core 
supernova in the sign at early times {t ~ 60 ms, around the 
neutronization peak): negative for ONeMg-core supernovae 
and positive for Fe-core ones. 
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If 013 is small, the Earth effect is non-vanishing at all 
times (fig. |7|, because the high-density resonance re- 
mains completely non-adiabatic. The double sign flip 
discussed above happen in this case as well, implying the 
same possibility for amplitude modulation and/or sign 
changes as for the scenario with intermediate 6*13. The 
major difference with respect to a Fe-core supernova is 
the different sign of the Earth effect at i ~ 60 ms, while 
at later times the two cases are similar. 
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FIG. 8: Same as fig. [ejfor sin^ 6*13 = 10"^. After the shock 
reaches the He shell, the conversion is essentially adiabatic 
and the Earth effect disappears. In a Fe supernova the effect 
would be zero at all times. 

For large ^13, oscillations in the Earth practically dis- 
appear with the shock passage through the high-density 
resonance (fig. |8|, due to the adiabatic character of this 
resonance in the post shock regime. This differs strongly 
from the case of a Fe-core supernova, where for the same 
values of ^13 (and the same normal mass hierarchy) the 
Earth effect is zero at all times. 

The sudden or progressive disappearance of the Earth 
oscillations is unique of a ONeMg-core supernova. For a 
Fe-core supernova the effect is either constantly non-zero 
or it appears at late times (i ^ 5 s) due to the shock 
wave propagation through the resonance layer [3]. 

We find similar results for other directions of propa- 
gation inside the Earth. The Earth effect is stronger for 
deeper trajectories inside the Earth, reaching about 40% 
size at 70 MeV for neutrinos moving along the diame- 
ter of the Earth. 



B. Spectra 

Finally, we find it useful to show the energy spectrum 
of the electron neutrino flux in a detector for a ONeMg- 
core supernova. It is given in the left column of fig. [TO] 
for different times and different values of 6113 . The right 
column of the figure illustrates the same spectra, but with 
the oscillation effects calculated for a Fe-core supernova 
(from E) mi- 



The figure shows two dramatic signatures of a ONeMg- 
core supernova in the high luminosity flux from neu- 
tronization {t ~ 60 ms), that were already pointed out 
in [7j : the lack of suppression of the flux, regardless of 
the value of 6*13, and the step in the spectrum at i? ~ 12 
MeV, due to neutrino-neutrino scattering. The exact po- 
sition and shape of the step in the energy spectrum would 
require a full detailed calculation that is not done here. 
However, using the results in [2S] we have checked that 
these quantities are rather independent of the details of 
the original neutrino spectra and of the matter density 
profile, and this assures the validity of fig. [TOjfor illus- 
tration. 

A third, well visible, feature is the negative sign of the 
Earth effect in the neutronization peak for a ONeMg-core 
supernova in contrast with the positive sign for a Fe-core 
one. This could be precious to distinguish the two pro- 
genitor types if 013 is small or unknown: indeed for small 
013 the neutronization peak is unsuppressed for both 
types and the step in the spectrum for the ONeMg-core 
supernova could be masked by the small statistics and 
poor energy resolution if it is near the threshold of the 
detector (e.g., a 7 MeV threshold for a water Cherenkov 
detector). The Earth effect, instead, is largest at high 
energy, where the energy resolution is better. 

The features discussed for a ONeMg-core supernova 
in the post-shock phase, namely the progressive decrease 
of the survival probability and of the Earth effect for 
intermediate 0i3, are present, but not visible at the scales 
used in fig. 10 They are of the order of ~ (5 — 25)% 
depending on the energy and of the specific value of ^13. 

Notice that, while p can decrease with time by a factor 
of 2 or more (fig. |4]), the net effect on the neutrino spec- 
trum is modest. This is due to the fact that p is small in 
value: the observed Ve flux is dominated by the originally 
produced u^^, Vt fluxes, and therefore it is impacted only 
at the level of tens of per cent by even major changes in 
the surviving component. 



V. DISCUSSION 

ONeMg-core collapses amount to 4% to 20% of all su- 
pernovae in the local universe [8] . If one of them happens 
in our galactic neighborhood, the observed conversion ef- 
fects will be a unique way to confirm the presence of the 
step in density at the base of the He shell and the faster 
shock propagation relative to a Fe-core supernova. The 
main oscillation signatures for both supernova types are 
summarized in fig. [9] In more detail, the presence of the 
step characteristic of a ONeMg-core supernova would be 
confirmed if: 

• 013 is known to be large and the neutronization 
peak does not disappear. Disappearance is pre- 
dicted for the smooth profile of a Fe-core supernova, 
where the neutrino propagation is completely adi- 
abatic. 
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• 6*13 is intermediate, small, or unknown, and data in- 
dicate a Vf, survival probability larger than ^ 0.32 
in the first 100-150 ms of the burst. This indi- 
cates that at least part of the produced z/g's exit 
the star in the state v\, which can not happen for 
the smooth density profile of a Fe-core supernova. 

• The Earth effect is negative in coincidence with the 
ncutronization peak (t ~ 60 ms). This is another 
sign that z/g's are converted into v\S. 

• Shock effects - regardless of when they happen ~ 
evidence a change from less adiabatic to more adi- 
abatic conversion as time passes. Such a change 
manifests itself as decrease of the i^e survival prob- 
ability (increase of Vf.-v^ permutation) and decrease 
in the amplitude of the oscillations inside the Earth. 
In absence of the step, the shock passage would 
have the opposite effect [T]. 

The timing of the shock effects will test the scenario 
of the faster shock propagation unambiguously. These 
effects are characterized by features that move through 
the neutrino spectrum, from low to high energy, as time 
passes. Such features can, in principle, provide a valuable 
measurement of the speed of the shock. 



of oscillations in the Earth drop quickly to zero in 
the post-shock phase (within ~ 300 ms). 

• An intermediate 6*13 (sin^ Bxz^ 3 • 10"^ - 3 • 10"^) is 
identified by the fact that both p and the amplitude 
of oscillations in the Earth decrease more slowly, 
and may vanish only in the late post-shock phase, 
t 700 ms. This late transition will begin with the 
lower energy neutrinos and extend to higher energy 
as time passes. The probability p changes in three 
distinct steps corresponding to pre-shock, early and 
late post-shock. For sin^ 6*13- 3 • 10"" - 3 • 10"^ a 
similar three-steps behavior characterizes the Earth 
oscillations at least in part of the energy spectrum 
(fig- §• 

• A small 6*13 (sin^ 6*13^ 3 • 10~^) will appear from 
the fact that p never vanishes, but only drops from 
p ~ 0.68 to p ^ 0.32 after the shock reaches the 
base of the He shell. The Earth oscillations also re- 
main present at all times, either unchanged or with 
subtle time modulations. These can be a tempo- 
rary flip of sign and/or a temporary reduction or 
even disappearance of the amplitude. They are due 
to the change of Pl from 1 to 0, passing by the crit- 
ical value of Pl = 1/2, for which the Earth effect 
disappears. 
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FIG. 9: Summary of the oscillation signatures of ONeMg-core 
supernovae compared to Fe-core supernovae at different times. 
They refer to the i/^ channel for the normal mass hierarchy; p 
is the Vs survival probability. The shaded areas represent the 
pre-shock phase, defined as the time interval before the shock 
front reaches the location of the MSW resonances, while the 
white regions ("shocked") indicate the later times. 

The time-dependent features induced by the shock 
wave on the observed signal will also give information 
on 013 . In particular: 

• A large 6*13 (sin^ ^13^ 3 • lO'^) will be singled out if 
both the survival probability p and the amplitude 



In the absence of Earth shielding, one would need to 
measure the value of the survival probability p and look 
for its time variations. Both could be difficult because 
the neutrino ffuxes in the different ffavors at the neutri- 
nosphere are uncertain and evolve with time, a fact that 
can mask time features due to ffavor conversion. The 
non-disappearance of the neutronization peak and step- 
like features that move from low to high energy, like that 
in Fig. |4] would probably be the easiest to distinguish. 

If a detector is shielded by the Earth, instead, the 
chances to distinguish the early shock effects are much 
better, since the oscillatory distortions of the neutrino 
spectrum induced by the Earth can not be mimicked by 
any other phenomenon |56l I57j . and their phase is well 
known thanks to the relatively precise measurements of 
the parameters Am^i and 9i2- 

If a supernova is optically obscured, neutrino data will 
be the only source of information on the nature of the 
progenitor. In such an event, searches for the oscilla- 
tion signatures of an ONeMg-core supernova would have 
an even higher importance. Information from data at 
early times will be complemented and substantiated by 
the data at late times, i ^ 5 s, where shock effects are 
expected for a Fe-core supernova and ^13 in the interme- 
diate and large range [21 IH [5l |6] . For example, the late 
appearance of oscillations in the Earth would exclude an 
ONeMg-core supernova. For small ^13 any shock wave 
signature disappears for a Fe-core supernova, making 
it more challenging to distinguish between Fe-core and 
ONeMg-core types. Discrimination is still possible from 
looking for signs of non-zero in the pre-shock phase. 
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which is typical of an ONeMg-core supernova. These are 
a large v^. survival probability, p ^ 0.32, and the negative 
Earth effect around the neutronization peak. 

In any event, the oscillations discussed here will be 
crucial for the correct interpretation of data from an 
ONeMg-core supernova, aimed at reconstructing the 
original neutrino fluxes, which are so important to test 
the theory of core collapse, reconstruct the mass of the 
neutron star and - indirectly - its equation of state, and 
to discuss the conditions for supernova nucleosynthesis. 

Before concluding, a word of caution is necessary about 
the validity of our results. They were obtained assum- 
ing that the conversion of neutrinos in the star can be 
described in terms of transitions between the eigenstates 
of the Hamiltonian - which justifies Eq. ^ - with the 



transition probabilities as discussed in sec. Ill While 
all the literature available at this moment supports this 
prescription, a confirmation from a detailed numerical 
study of the problem still lacks. We expect that, after 
such study is done, our results remain valid at least in 
the main message that a distinctive pattern of neutrino 
conversion is associated with a ONeMg-core supernova, 
characterized by high survival probability of v,, due to the 
step in density at the base of the He shell and by early 
shock effects that drive the conversion towards stronger 
fiavor conversion. 

Our work is focused exclusively on the case of neutri- 
nos and normal mass hierarchy, for two main reasons. 
The first is that the case of normal hierarchy is particu- 
larly simple in its oscillation pattern, and has also been 
studied in the most detail for ONeMg-core supernovae in 
the pre-shock phase [7] . This makes it an ideal choice for 
an initial study of the oscillation signatures of a ONeMg- 
core supernova. The second motivation is that the neu- 
trino channel is especially interesting for a ONeMg-core 
supernova because of the very distinctive signature in 
the neutronization peak, i ~ 60 ms (fig. 10). General 
arguments on neutrino conversion in the dense matter 
of a supernova (see e.g. [Uj) suggest that the features 
expected for neutrinos and normal hierarchy should ap- 
pear - same in character but different in the details ~ 
also in the other cases of antineutrinos and/or inverted 
mass hierarchy. Specifically, for the inverted hierarchy 
all the phenomena associated with the high density res- 
onance will appear in the antineutrino channel. Also, 
neutrino self-interaction will induce strong fiavor conver- 
sion on both neutrinos and antineutrinos within the first 
few hundreds of Km of radius (see e.g. HO]). 

The generalization of our results to other combina- 
tions of channels (neutrinos and antineutrinos) and hier- 
archies require a more in-depth modeling of the effects of 
neutrino-neutrino scattering, a difficult task that receives 



further motivation from this work. The combination of 
observations in different channels can only strengthen the 
conclusions that a single channel can give. 

This brings us to the question of what detector is best 
for the study of the effects we have discussed. For the case 
considered here (neutrino channel), the optimal setup is 
an experiment that can detect electron neutrinos exclu- 
sively. Liquid argon technology is very suitable for this. 
To have sufficient statistics in time bins ~ 100 ms wide, 
necessary to test the fast time variations found here, large 
volumes are required. Thus, the massive liquid argon 
project MODULAr (20 Kt mass, corresponding to ^ 10"^ 
Ve events from a galactic supernova [59^) and the even 
more ambitious GLACIER ^6D] (up to 100 Kt mass, i.e. 
~ 6 • 10^ Ve events) would be very valuable. Consider- 
ing that no large Vf, detector is currently active, after the 
closure of SNO [6lj , it is extremely important to empha- 
size the importance of having new detectors of this type 
to study supernovae, and our paper contributes in this 
direction. Of course, the best out of a supernova ob- 
servation would be obtained from combining data from 
different detection channels. The effect on antineutri- 
nos due to the step in the density of an ONeMg-core 
supernova would be very well visible in the inverse beta 
decay events at SuperKamiokande (see e.g., [62]) and at 
its larger versions, UNO, HyperKamiokande and MEM- 
PHYS [Ml [Ml , as well as at the planned liquid scin- 
tillator experiment LENA [66^. With several detectors 
running, there is a substantial chance that at least one of 
them will be shielded by the Earth I5H HT] and thus will 
have a enhanced sensitivity to neutrino flavor conversion. 
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FIG. 10: Left column: predicted spectra of the i/g flux in a detector for ONeMg-core supernova, inclusive of oscillations in 
the star and in the Earth (nadir angle 60 degrees). Right column: the same spectra, but with the oscillation effects that are 
characteristic of a Fe-core supernova. From the upper to lower panel: sin'^ 613 = 0.01,6 • 10~'',10~^. The thick curves refer 
to different times: t = 60,450,700 ms (solid, short dashed, and long dashed, respectively). For t = 60 ms we also show the 
spectrum in absence of Earth shielding (thin solid line). The vertical axis has units of MeV~^s~^. 



